Sequential and simultaneous emission of fragments from p+Al collisions 
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The energy and angular dependence of double differential cross sections d 2 a/dQdE were measured 
for p, d, i, 3 ' 4 ' 6 He, 6 ' 7 ' 8 Li, 7 ' 9 ' 10 Be, and 10 ' u ' la B produced in collisions of 1.2, 1.9, and 2.5 GeV protons 
with an Al target. The spectra and angular distributions of Li, Be, and B isotopes indicate a presence 
of two contributions: an isotropic, low energy one which is attributed to the evaporation of particles 
from excited remnants of the intranuclear cascade, and an anisotropic part which is interpreted to 
be due to multifragmentation of the remnants. It was found that such a hypothesis leads to a very 
good description of spectra and angular distributions of all intermediate mass fragments ( 6 He - 12 B) 
using the critical value of the excitation energy per nucleon as a single parameter, varying slowly 
with the beam energy. 
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I. INTRODUCTION 

In recent studies of reactions induced by GeV protons 
on Au [J [5] and Ni targets [21 0], it has been found 
that the inclusive spectra of intermediate mass fragments 
(IMF), i.e., particles with Z> 3, but lighter than fission 
fragments, as e.g. Li, Be, B, etc., contain two compo- 
nents which differ in energy and angular dependencies. 
The low energy component of the spectra is almost angle 
independent, while the high energy part of the spectra 
changes with the angle becoming more steep with in- 
creasing scattering angle. These properties of the spec- 
tra could not be quantitatively reproduced by the tradi- 
tional two step model which assumes that the first stage 
of the reaction proceeds via an intranuclear cascade of 
the nucleon- nucleon and nucleon-pion collisions leaving 
the excited target remnant in equilibrium whereas the 
second stage consists in the evaporation of nucleons and 
composite particles. 

The same effect of the presence of two components in 
the spectra has been observed for various target nuclei by 
other authors, e.g., by Green et al. [51 [5] for reactions in 
p+Ag system at 2^=0.21, 0.3, and 0.48 GeV, by Herbach 
et al. 7 for target nuclei between Al and Th at 1.2 
GeV, by Letourneau et al. [8] for p+Au collisions at 2.5 
GeV, by Westfall et al. for C, Al, Ag and U targets 
irradiated by protons of 2.1 and 4.9 GeV energies, by 
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Hyde et al. [TO] for p+Ag system at 5.5 GeV, and by 
Poskanzer et al. [TJj for p+U reactions at 5.5 GeV. 

The analysis of the energy and angular dependencies of 
differential cross sections for IMFs from p+Au and p+Ni 
reactions PHI] has shown that the data can be well repro- 
duced by a phenomenological model which assumes that 
the composite particles are emitted isotropically from two 
sources moving along the beam direction. The break-up 
of bombarded nuclei, occuring in the first stage of the 
collision, was proposed [JJ to account for the observed 
phenomenon. 

The competitive hypothesis, which is studied in the 
present work, states that the excited nuclei are divided 
into two groups not due to the geometrical conditions 
imposed by their creation mode but on the contrary, by 
their decay. It may be conjectured that the decay mode 
of the excited nuclei depends mainly on their excitation 
energy per nucleon. The nuclei which have small excita- 
tion energy evaporate nucleons and composite particles, 
whereas the highly excited nuclei (above some critical 
energy per nucleon) may undergo a phase transition, i.e. 
multifragmentation appears [321 H3]. The energy avail- 
able in the multifragmentation is usually larger than that 
in the evaporation process because the multifragmenta- 
tion appears only for highly excited nuclei. Furthermore, 
the full excitation energy is then released in a single act, 
whereas the evaporation consists of several stpdf in which 
light particles take a part of the excitation energy and 
heavy residua of the evaporation have smaller energy - 
due to the momentum conservation in each consecutive 
two-body decay of the de-exciting nucleus. Thus, the 
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intermediate mass fragments emitted from a multifrag- 
mentation could have on average larger energies than the 
evaporation residua with the same mass. Due to this 
fact, the evaporation should give rise to qualitatively dif- 
ferent spectra than emission of the same heavy particles 
in the multifragmentation. This effect can explain the 
appearance of two qualitatively different components of 
the spectra of intermediate mass fragments. Thus, the 
present study concentrates on the investigation of the 
discussed above hypothesis concerning properties of the 
spectra of intermediate mass fragments. 

Since exact calculations of multifragmentation of heavy 
nuclei are practically not possible due to an enormous 
number of involved partitions, and since the intranuclear 
cascade models, which neglect the shell structure of nu- 
clei might not work well for very light nuclei (cf., e.g., 
Z. Fraenkel et al. [2]), the choice of the target nucleus 
for investigation of validity of the above proposed hy- 
pothesis should fulfill a compromise between these two 
constraints. The aluminium target has been therefore 
selected for the present study. 

The paper is organized as follows: the experimental 
details are discussed in the next section, the analysis of 
the IMF data performed by means of the intranuclear 
cascade (for the first stage of the reaction) combined with 
evaporation and multifragmentation (for the second step) 
is presented in the third section, and the last section 
summarizes the obtained results. 



II. EXPERIMENTAL RESULTS 

The PISA experiment has been performed using the 
internal beam of COSY (COoler SYnchrotron) of the 
Research Center in Juelich. The apparatus and exper- 
imental procedure have been described in previous pub- 
lications [TJ [H 0] thus in the present work only details, 
characteristic for the studied reactions are discussed. 

Self-supporting aluminium target of the 170 /ig/cm 2 
thickness was bombarded by the internal proton beam 
of COSY. Three beam energies were used: 1.2, 1.9, and 
2.5 GeV. To assure the same experimental conditions for 
all beam energies COSY operated in the so called super- 
cycle mode. In this mode several cycles were alternated 
for each requested beam energy, consisting of protons in- 
jection from the cyclotron JULIC to COSY ring, their 
acceleration with the beam circulating in the ring below 
the target, and irradiating the target by slow movement 
of the beam in the upward direction. Due to the appli- 
cation of the supercycle mode of the target irradiation 
all conditions of the experiment except the energy of the 
proton beam remained unchanged. This allowed to min- 
imize the effect of systematic uncertainties on the energy 
dependence of the measured cross sections. 

Double differential cross sections d 2 <j/d£ldE were mea- 
sured at seven scattering angles: 15.6°, 20°, 35°, 50°, 
65°, 80°, and 100°. The mass and charge identification 
of detected particles was realized by the A E - E method 



using telescopes consisted of silicon semiconductor detec- 
tors backed (in four cases: 15.6°, 20°, 65°, and 100°) by 
a 7 cm thick Csl detector with a photo-diode readout, 
which were used to detect high energy light charged par- 
ticles passing through the silicon detectors. 

The d 2 a/dttdE were measured for the following ejec- 
tiles: p, d, t, 3 ' 4 < 6 He, 6 < 7 ' 8 < 9 Li, 7 < 9 > 10 Be, ^,U,12 B and c 
Here we will concentrate on the IMF data. Discussion 
on the light charged particle emission will be published 
separately. 

In order to allow for absolute normalization of differen- 
tial cross sections we performed a "two moving sources" 
fit to d 2 ajd£ldE for 7 Be production in our experiment 
and compared the resulting total cross section with data 
from ref . [TS] . It was found that the angular and energy 
dependence of d 2 a/d£ldE can be well reproduced by a 
simple formula representing the isotropic emission from 
two sources moving forward along the beam direction. 
Each source emitting particles with Maxwellian energy 
distribution was characterized by its velocity /3, temper- 
ature T, height of the Coulomb barrier between 7 Be and 
the source remnant described by parameter k, and by 
emission intensity a (see Appendix of ref. [T] for de- 
tails of the parametrization). The a parameter has the 
meaning of the energy and angle integrated cross sec- 
tion attributed to a given source. Thus, the not normal- 
ized cross section a a . u . for 7 Be production is equal to the 
sum of parameters u\ + a 2 for both sources. Best values 
of the parameters were found by fitting simultaneously 
the full set of the 7 Be spectra (seven scattering angles). 
Alas, the fits lead to ambiguous results because the ex- 
perimental spectra did not cover the full energy range 
allowed by kinematics. Low energy particles were not 
registered because of finite low-energy detection thresh- 
olds of the telescopes. This lack of information on the low 
energy part of spectra could strongly influence the value 
of the energy integrated cross section since the spectra 
have Maxwellian shape with its maximum lying in the 
neighborhood of the energy threshold. Fortunately, it 
turned out that the spread of values of a parameters was 
smaller than 10 % among the sets of parameters which 
provided the same, best % 2 values obtained for various 
combinations of fixed Coulomb barrier parameters k\ and 
&2 • The final values of the a parameters for both sources 
were taken as the arithmetic mean of results obtained for 
equivalent quality fits, i.e., those which have the same, 
smallest % 2 value. The error of the normalization factor 
obtained in such a way was w 9% for all studied energies. 
This value does not include the inaccuracy of the liter- 
ature value of the production cross section er( 7 Be) [T5] 
which is believed to be smaller than 10%. 

The good quality of the absolute normalization of the 
present experiment is confirmed by a comparison of the 
obtained differential cross sections with the data from 
literature. To our knowledge only one measurement of 
light charged particle spectra as well as intermediate 
mass fragment spectra is present in the literature for alu- 
minium target at proton beam energies similar to those 
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FIG. 1: Spectra of 4 He, 7 Li, 9 Be, n B measured at 35° for 
p+Al collisions at three proton beam energies: 1.2 (circles), 
1.9 (full circles), and 2.5 GeV (triangles). To avoid over- 
lapping of symbols only different energy bins are shown for 
different beam energies. 



used in the present study. This is the paper of Westfall 
et al. 9J dealing with reactions in the p+Al system in- 
vestigated at 4.9 GeV proton energy. This energy is 2-4 
times higher than those used in the present work, how- 
ever, it is known that the total production cross sections 
for this target vary only slightly at proton beam energies 
larger than 1 GeV (cf. e.g. ref. [T5]). The differential 
cross sections measured in the present experiment also 
change only slightly with the beam energy as it is shown 
in Fig. [T] As can be seen, the cross sections increase very 
slightly with the increasing beam energy. This increase is 
stronger for lighter ejectiles ( 4 He, 7 Li), whereas the data 
for heavier IMFs are in the limits of errors independent 
of the beam energy. 

The weak energy dependence of the differential cross 
sections leads to the conclusion, that it is reasonable to 
compare the present data with results obtained by West- 
fall et al. [5] at even higher incident proton energy. Such 
a comparison is depicted in Fig. [2] for lithium, beryllium 
and boron isotopes. As can be seen, the shapes and mag- 
nitudes of the spectra are in excellent agreement for all 
ejectiles. 
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FIG. 2: Experimental spectra of Li, Be, and B from the 
present work (circles) measured at 100° in p+Al collisions at 
proton beam energy 2.5 GeV and data from ref. [5] measured 
at 90° at proton energy 4.9 GeV (lines). 



III. THEORETICAL ANALYSIS 

The model calculations of the first stage of the reaction 
were done by means of the INCL4.3 computer program 
|16[ I17j realizing the intranuclear cascade of nucleon- 
nucleon collisions. The decay of remnant nuclei (A,Z) 
excited to the energy E* was evaluated in the frame of 
two different models: the generalized evaporation model 
GEM2 of S. Furihata [HlQlj and the multifragmentation 
treated as the Fermi breakup model using the computer 
program ROZPAD of A. Magiera [2"D] . 

The Fermi breakup model [2"lTf2"3"] assumes that the 
probability W of disintegration of an excited nucleus with 
mass M and volume V into channel containing n frag- 
ments with masses m, and spins Sj (i — 1,2, ...,n) is 
proportional to its phase space volume [23] 
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where g is the number of spin projections of n fragments 
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g = Y[ (2s-i + 1), G = 1/ Yi n r accounts for identity of 

i=l j=l 

k 

rij fragments of kind j (fc is defined by n = Y) rij), and 

i=i 

r(x) is the gamma function. The total kinetic energy of 
n fragments at the moment of breakup Eki n is related to 
the excitation energy E* of the decaying nucleus (A,Z) 
according to the equation: 



= {E* + Mc 2 ) - 



i=l 



rriiC 2 - U° 



where [/„ represents the Coulomb interaction energy of 
the fragments. 

The double differential cross sections cPa/dfldE are 
obtained in the Fermi breakup model by random gener- 
ation of energies and flight directions of fragments over 
the whole accessible phase space (23j . 

It should be pointed out that the Fermi breakup model 
gives statistical method of dealing with the multifrag- 
mentation process based on the assumption of full equi- 
libration of the excited system in which the energy and 
linear momentum conservation is exactly taken into ac- 
count. In the following the "multifragmentation" will be 
understood as such a specific model of this process. 

It was assumed that the sequential emission of IMFs, 
i.e., the evaporation of a single particle in each step, dom- 
inates when the remnants of the cascade are excited to 
low energies whereas the remnants with higher excita- 
tion energies split simultaneously into many pieces thus 
they are subject of multifragmentation. The critical ex- 
citation energy per nucleon (E* /A) cr , i.e., the smallest 
value of excitation energy per nucleon at which multifrag- 
mentation appears was treated in the present work as a 
free parameter. The only free parameter of the Fermi 
breakup model - the so called freeze- out radius rg was 
also fitted. It has the meaning of a reduced radius of 
the decaying, spherical nucleus of the mass number A, 
i.e. the volume V of the decaying nucleus is calculated 
as V — 47rrQ^4/3. The default value of the ro parameter 
used in the GEANT4 is equal to 1.4 fm [23]. 

The free parameters of the present model: freeze-out 
radius ro and critical excitation energy (E*/A) cr were 
searched for by comparison of theoretical spectra with 
the experimental data, for 6 Hc, 6 ' 7 > 8,9 Li, 7 > 9 ' 10 Be, and 
10 ' n ' 12 B. Chi-square values x 2 were evaluated summing 
over all these data, i.e. over all ejectiles listed above and 
over seven scattering angles for each of these ejectiles for 
several fixed values of ro parameter and several values of 
critical excitation energy (E*/A) cr . The obtained values 
of the chi-square are shown in Fig. [3] for proton beam 
energy 1.2 GeV. 

The following properties of the chi-square dependence 
are clearly visible: 

• Broad minima of the chi-square treated as a func- 
tion of the critical excitation energy (E* /A) cr are 
present for each fixed value of the r parame- 
ter. The chi-square increases very strongly when 
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FIG. 3: Results of the analysis for the proton beam energy 
equal to 1.2 GeV. Upper panel: chi-square values versus crit- 
ical value of the excitation energy per nucleon (E* /A) cr for 
several values of reduced freeze-out radius ro (in fm) . Lower 
panel: Minimal values of the chi-square for each value of ro. 



the critical excitation energy decreases to values 
smaller than « 5 MeV/nucleon. This behavior 
points to the fact that the nuclei at such low excita- 
tion energies are not subject of the multifragmen- 
tation but rather sequentially evaporate particles. 

• Minimal value of the chi-square for each value of 
the r parameter, found from the chi-square depen- 
dence on (E* /A) cr , is presented as a function of ro 
in the lower panel of the Fig. [3] . As can be seen, 
this dependence may be well approximated by a 
concave parabolic function, which minimum allows 
to choose the best fit value of the ro parameter. 

The fit procedure described above has been applied to 
the data measured at all three proton beam energies: 1.2, 
1.9, and 2.5 GeV, leading to the same (in the limits of 
errors) best fit values for the reduced radius ro: 1.40(31) 
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FIG. 4: Experimental spectra (circles) measured at three scat- 
tering angles: 20°, 50°, and 100° (left, middle, and right 
columns) for 6 Li, 7 Li, 8 Li, and 9 Li (upper, middle, and lower 
rows) for p+Al collisions at proton beam energy of 1.2 GeV. 
The lines represent theoretical calculations: blue (the steep- 
est), red (the least steep), and black lines correspond to evap- 
oration, multifragmentation, and their sum, respectively. 



fm, 1.40(25) fm, and 1.39(43) fm for beam energy 1.2, 
1.9, and 2.5 GeV, respectively. It is worth pointing out 
that the obtained values of rg parameter are equal in 
the limits of errors to tq =1.4 fm assumed as the default 
value in the literature |23j . Thus fixing the parameter 
ro at ro = 1.40(20) fm allows to find the best fit value 
of the second parameter - the critical excitation energy 
(E* /A) cr - independently for each beam energy. 

It was found that the critical excitation energy is also 
almost the same in the studied beam energy range. It 
is equal to 7.7(8) MeV/nucleon, 7.0(8) MeV/ nucleon, 
and 6.1(9) MeV/nucleon for 1.2, 1.9, and 2.5 GeV beam 
energy, respectively. 

The quality of the description of experimental data for 
p+Al collisions may be judged from inspection of Figs [4] 
- [6j where a representative sample of the data measured 
at proton beam energy equal to 1.2 GeV is shown. As 
can be seen, the contribution of multifragmentation to 
the spectra (shown by the least-steep, red line in the 
figures) is crucial for reproducing high energy parts of 
the spectra of IMFs. On the other hand the evaporation 
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FIG. 5: Same as in Fig. [Zjbut for 7 Be, 9 Be, and 10 Be, respec- 
tively 

which leads to very steep shape of the spectra is not able 
to describe their high energy tails but dominates at low 
energies of the ejectiles. It turned out that the sum of 
evaporation and multifragmentation describes very well 
the full spectra of 6 ^ 7 > 8 ' 9 Li, 7 - 9 ' 10 Be, and 10 > n > 12 B.' Note 
that the fluctuations of the theoretical curves visible in 
the figures have no meaning - they are only due to a 
limited statistics of the Monte Carlo calculations. 

It should be emphasized that all the spectra shown 
here, obtained at 1.2 GeV proton beam energy, are re- 
produced assuming the same values of the parameters: 
r = 1.4 fm and (E*/A) cr = 7.7 MeV. 

An equally good description of IMF data was achieved 
for higher beam energies (1.9 and 2.5 GeV) using the 
same value of r$ parameter and appropriate critical ex- 
citation energies. Since the shape of the spectra almost 
does not change with the beam energy, the Figs [4] - [6] 
are representative for the quality of data reproduction 
obtained at all studied energies. 

The angle and energy integrated double differential 
cross sections cPa/dfldE obtained in the analysis de- 
scribed above are collected in Table |TJ The following con- 
clusions can be derived: 

• Evaporation cross sections decrease quickly with 
beam energy for all IMFs, on the contrary to 
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FIG. 6: Same as in Fig. [Zjbut for 10 B, n B, and 12 B, respec- 
tively 

increasing of multifragmentation cross sections. 
These monotonic dependencies lead to the fast in- 
crease of the relative contribution of multifragmen- 
tation to the total cross section. 

• The multifragmentation contribution dominates for 
all ejectiles at the highest beam energy. It exhausts 
at least 51% of the total cross section (for n B) but 
for a few ejectiles it reaches even 95% (for 8 Li and 
9 Li). 

• The sum of both contributions increases with the 
beam energy but not as quickly as the fragmenta- 
tion cross sections themselves because of decreasing 
of the evaporation contribution. The ratio of the 
total production cross section at the beam energy 
of 2.5 GeV to that at 1.2 GeV varies from about 
1.1 for 6 Li and n B to about 2.2 for 9 Li. 

It was checked whether the parameter of the present 
model, i.e. the critical excitation energy per nucleon 
(E* /A) cr is compatible with the critical excitation en- 
ergy found from the analysis of the caloric curves. For 
this purpose the values of the critical energy obtained in 
the present investigation were compared with those from 
a compilation of caloric curve results published by Na- 
tovitz et al. [23]. The target mass dependence of the 



TABLE I: Total production cross sections of intermediate 
mass fragments for p+Al collisions at three proton beam en- 
ergies: 1.2, 1.9, and 2.5 GeV. In the first column the symbol 
of the ejectile is listed, in the second column the beam en- 
ergy is specified, and in the following columns the production 
cross section due to the sequential evaporation, production 
cross section due to the multifragmentation, the sum of both 
cross sections, and the relative contribution of the multifrag- 
mentation to the total cross section are presented. 
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critical excitation energy obtained from the analysis of 
the caloric curves by Natovitz et al. is presented in Fig. 
[7] together with critical energies found from p+Al colli- 
sions investigated in the present work. The former results 
were averaged in Ref. [7] over thirty units of mass of de- 
caying nuclei and were attributed to the central value 
of each mass region. The present results correspond to 
multifragmentation of an ensemble of excited nuclei with 
the average mass < A >= 18.7 and with the standard 
deviation of masses <j(A) = 2.6. It is clear that the 
present values of the critical excitation energy fit per- 
fectly to the compilation of data obtained from the study 
of caloric curves. The values of the critical energies from 
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FIG. 7: The full squares represent the compilation of critical 
excitation energies per nucleon as a function of mass number 
of decaying nuclei obtained from the study of caloric curves 
|24| , the solid straight line shows the linear regression of these 
data whereas the dashed and solid hyperbolas correspond to 
one- and two- standard deviation confidence intervals of the 
regression line, respectively. The empty symbols depict the 
critical excitation energy values used in the present analysis 
for p+Al collisions; circle, square, and triangle correspond to 
1.2, 1.9, and 2.5 GeV beam energy, respectively. 



the present study lay for all three proton beam energies 
inside the confidence interval of one standard deviation 
around the regression line approximating the mass de- 
pendence of the critical energy per nucleon. 



IV. SUMMARY 

In the present work the hypothesis was investigated, 
claiming that the energy and angular distributions of the 
double differential cross sections d 2 ajd£ldE for interme- 
diate mass fragments from inclusive measurements of 
reactions proceeding in the p+Al system at GeV proton 
beam energies can be reproduced by a reaction model in 
which a specific two-step mechanism is involved. In the 
first stage of the reaction a cascade of nucleon-nucleon 
and pion-nucleon collisions appears leaving an ensemble 
of the excited residual nuclei. These nuclei decay in the 
next stage of the reaction according to two competing 
scenarios: (i) The nuclei with the excitation energy per 
nucleon smaller than some critical value (E*/A) cr - 
treated as a parameter, decay by a sequential emission 
of evaporated particles, (ii) the nuclei excited to energies 
above the limiting value (E*/A) cr are subjects to 
multifragmcntation. The first of these mechanisms is 
responsible for populating the low energy part of the 
spectra, whereas the second leads to emission of also 
higher energy intermediate mass fragments. 



The excellent agreement of the calculated within such 
a model differential cross sections d 2 a/dildE with the 
experimental data confirms the validity of the assumed 
hypothesis. 

The investigations performed in the present work pro- 
vide also the following conclusions: 

• It was found that all intermediate mass fragments 
are products of the second stage of the reaction 
with a relative contribution of multifragmentation 
of the remnants of the nucleon-nucleon cascade 
systematically increasing with the proton beam 
energy. 

• The only free parameter of the calculations men- 
tioned above was the critical excitation energy 
(E* I A) cr which decides which option of the decay 
of the excited remnants of the nucleon-nucleon 
cascade should be chosen. The second parameter 
which could in principle vary, i.e. the reduced 
freeze-out radius ro, appeared to be independent 
of the beam energy and to be equal to 1.4 fm, i.e. 
to the default value used in the Fermi breakup 
model (see e.g. ref. |23|). 

• The value of the critical energy found from the 
present analysis of the inclusive spectra is com- 
patible with those obtained from the study of the 
caloric curves. Therefore, it may be concluded that 
the investigation of intermediate mass fragment 
spectra obtained in inclusive measurements may 
be used as an alternative method to studying of 
the caloric curves for the extraction of the critical 
excitation energy per nucleon. 

Successful description of the intermediate mass frag- 
ment emission from p+Al collisions in the 1.2 - 2.5 GeV 
proton beam energy range suggests that a similar picture 
of the reaction mechanism might be also appropriate for 
heavier nuclear systems as, e.g., p+Ni or p+Au where 
another model of the IMF emission has been proposed 
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